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Abstract 

Genetic differentiation of morphologically variable members of 
the G. pulex- group in northwestern Europe was investigated by 
electrophoresis at 20 enzyme loci. Five recently described related 
forms were examined with reference to the “classical” species 
G. pulex , G. wautieri and G.fossarum. G. pulex and G. wautieri 
were shown to be genetically homogeneous and clearly distinct 
species. In the other species, morphologically similar forms were 
shown to be genetically distant, and inter-areal intraspecific 
genetic differentiation may warrant recognition of subspecies or 
sibling species. In some cases, levels of reproductive isolation of 
the forms involved could be assessed by cross-breeding experi¬ 
ments. To estimate dispersal capabilities, gene flow levels were 
indirectly estimated from gene frequency data. Gene flow levels 
are generally low and fluctuating population sizes may enhance 
the occurrence of stochastic processes. Dendrograms derived 
from genetic distances were compared with zoogeographical and 
paleoclimatological evidence. The distribution areas of the spe¬ 
cies involved correspond to patterns predicted by these data. 

Resume 

La differentiation genetique des membres morphologiquement 
variables du groupe Gammarus pulex dans le nord de l’Europe 
occidentale a 6t6 etudiee par electrophorese k 20 loci d’enzymes. 
Cinq formes apparentees recemment d6crites ont 6te etudiee par 
rapport aux especes “classiques” G. pulex , G. wautieri et G. 
fossarum. G. pulex et G. wautieri se sont months des espfcces 
genetiquement homogenes et clairement distinctes. Dans les 
autres especes, des formes morphologiquement similaires se sont 
montrSes genetiquement eloign^es, et la differentiation g£n£- 
tique intraspetifique entre aires de repartition diffSrentes peut 
permettre de reconnaitre des sous-espfcces ou des especes- 
jumelles. Dans certains cas, des essais d’hybridation ont permis 
d’etablir des niveaux d’isolation reproductive des formes con¬ 
cernees. Pour une estimation des capacity de dispersion, des 


niveaux de flux genetique ont ete indirectement estimes k partir 
des donnees de frequence des genes. En general, les niveaux de 
flux genetique obtenus sont bas et il est probable que Foccur- 
rence des processus stochastiques est favorisee par les fluctua¬ 
tions de Fordre de grandeur des populations concernees. Des 
dendrogrammes derives des distances genetiques ont ete con¬ 
fronts aux donnees zoogeographiques et paleoclimatologiques. 
Les aires de repartition des especes concernees correspondent 
bien avec les modeles ptdits par ces donnees. 

Introduction 

Due to a considerable amount of morphological 
variability, the systematic position of the species of 
the Gammarus pulex- group (sensu Karaman & 
Pinkster, 1977) in western Europe has been a source 
of confusion for quite some time (see Pinkster, 
1972). After careful morphological examination 
and hybridization experiments by several authors 
(e.g. Wautier & Roux, 1959; Roux, 1971; Goed- 
makers, 1972; Pinkster, 1972; Karaman & Pinkster, 
1977) most of the “deviating” forms encountered 
could be attributed to the “classical” species G . pu¬ 
lex (Linnaeus, 1758), G.fossarum Koch, 1835, and 
G. wautieri A.L. Roux, 1967. 

However, three forms were found that are closely 
related to, but clearly distinguishable from G. fos¬ 
sarum (cf. Goedmakers, 1974; Goedmakers & 
Roux, 1975; Pinkster, 1983; Pinkster & Scholl, 
1984; Scheepmaker, 1987). These three forms show 
morphological affinity with G. ibericus Margalef, 
1951, but are reproductively isolated from it (Pink¬ 
ster, 1983; Pinkster & Scholl, 1984). Two of these 
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Table /. Sampling localities and species distribution. 


Station 

no. 

Area in 
fig- 1A 

Species 

Country 

Prov/dept 

Drainage 

system 

Locality description/ reference 

Sampling 

date 

1 

F 

G. orinos 

France 

Ardeche 

Ardeche 

Rivulet draining into the Ardeche, at the 
crossing with road D243, 3.5 km S.E. of 
Chirols 

10-VII-86 

2 

F 

G. orinos 

France 

Ardeche 

Ardeche 

Brooklet draining into the Ardfcche, at 
crossing with road D536,0.1 km N. of 
Montpezat-sous-Bauzon 

10-VII-86 

3 

D 

G. orinos 

France 

Confeze 

Dordogne 

Brooklet draining in the Triouzoune, at 
crossing with road D30,2 km E. of Meymac 

13-VII-’86 

4 

D 

G. orinos 

France 

Contze 

Vienne 

Rivulet at the source of the Vienne, at crossing 
with road D21, 8 km S.E. of Peyrelevade 

13-VIK86 

5 

D 

G. orinos 

France 

Coirfcze 

Dordogne 

Ozanges, at the crossing with road D979, 

1.5 km N.E of St. Victour 

13-VIK86 

6 

E 

G. orinos 

France 

Carnal 

Dordogne 

Rivulet draining into the Maronne, alongside 
the railroad in St. Christophe-les-Gorges, 
near road D6 

15-VII-86 

7 

G 

G . cf. fossarum 

France 

Pyr6n6es 

Orientales 

Agly 

Rivulet draining into the Desix, at crossing 
with road D2,2 km W. of Soumia 

19-VII-'86 

8 

G 

G. cf. fossarum 

Fiance 

Pyr6n6es 

Orientales 

Agiy 

Rivulet draining into the Desix, at crossing 
with road D2,2 km W. of Rabouillet 

19-VIK86 

9 

C 

G. fossarum s.l. 

France 

Is£re 

Isfcre 

Vence, alongside road D57,0.5 km E. of 
Sarcenas 

8-VII-'86 

10 

C 

G. fossarum s.l. 

France 

Isfcre 

Is£re 

Mountain stream draining into the Is£re, at 
crossing with road D29, 0.5 km S.W. of 

St. Pancrasse 

8-VII-'86 

11 

C 

G. fossarum s.l. 
G. wautieri 

France 

Is&re 

Rhone 

Mountain stream draining into the Guiers, at 
crossing with road D520, 2 km W.N.W. of 

St. Pierre-de-Chartreuse 

7-VIK86 

12 

C 

G. fossarum s.l. 
G. wautieri 

France 

Is&re 

Is£re 

Mountain stream draining into the Isfcre, at 
crossing with road D29, 0.5 km S.W. of 

St. Bernard 

8-VII-86 

13 

B 

G. fossarum 

G. wautieri 

France 

Isfcre 

Rhone 

Amby, at the Optevoz laundry place, alongside 
roadD52A 

6-VII-86 

14 

C 

G. fossarum s.l. 
G. wautieri 

France 

Isere 

Is&re 

Tenaison, at the crossing of roads D105 and 
D105A, 2 km N. of Proveysieux 

9-VII-'86 

15 

F 

G. fossarum s.l. 

France 

Ard&che 

Ardeche 

Brooklet draining into the Ardeche, at crossing 
with road D323,2 km S.W. of Prades-les- 
Coutours 

ll-VII-'86 

16 

F 

G. fossarum s.l. 

France 

Ard&che 

Ardeche 

Brooklet draining into the Ardeche, at crossing 
with road D323,0.1 km W. of St Cirgues- 
de-Prades 

ll-VII-’86 

17 

■ 

G. fossarum s.l. 

Ranee 

Bouches-du- 

Rhone 

Rhone 

Sorgue, alongside road D25, in Fontaine-de- 
Vaucluse 

6-XI-*85 

18 

I 

G. fossarum s.l. 

France 

Bouches-du- 

Rhone 

Arc 

Station 13 in Scheepmaker (1987) 

29-VII-'86 

19 

I 

G. stupendus 

France 

Var 

Argens 

Station 17 in Scheepmaker (1987) 

22-VII-’87 

20 

I 

G. stupendus 

France 

Var 

Huveaune 

Station 43 in Scheepmaker (1987) 

29-VII-86 

21 

C 

G. wautieri 

France 

I sere 

Is£re 

Mountain stream draining into the Isfcre, 
alongside road D512,1 km S. of Le Sappey 
en Chartreuse 

8-VII- , 86 

22 

G 

G. p. gallicus 

France 

Pyr6n£es- 

Orientales 

Tet 

Spring at 0.1 km E. of road D46, 1.5 km 

S.S.W. of Millas 

6-VUI-'87 

23 

G 

G. p. gallicus 

France 

Pyr6n6es- 

Qrientales 

let 

Spring/brooklet crossing track alongside 
the river Tet, near road N116, 0.5 km E. of 
Nefiach 

is-viii-^ 

24 

I 

G. p. gallicus 

France 

Bouches-du- 

Rhone 

Arc 

station 16 in Scheepmaker (1987) 

5-XI-'85 
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Table 1. Cont. 


Station 

Area in 

Species 

Country 

Prov/dept 

Drainage 

Locality description/ reference 

Sampling 

no. 

fig. 1A 




system 


date 

25 

- 

G. p. pulex 

France 

Loiret 

Loire 

Clfry, at crossing with road D32,1.5 km 

N.W. of St. Hilaire-les-AndrSsis 

16-Vni-’88 

a 

F 

G. orinos 

France 

Ardfcche 

Ard&che 

Station 2 in Scheepmaker et al. (1988) 


b 

H 

G. cf. fossarum 

France 

Arifcge 

Arifcge 

Station 3 in Scheepmaker et al. (1988) 


c 

H 

G. cf. fossarum 

France 

Arifcge 

Ariege 

Station 4 in Scheepmaker et al. (1988) 


d 

H 

G. cf. fossarum 

France 

Arifege 

Garonne 

Station 5 in Scheepmaker et al. (1988) 


e 

I 

G. stupendus 

France 

Var 

Gapeau 

Station 1 in Scheepmaker (1987) and 

Scheepmaker et al. (1988) 


f 

M 

G. fossarum s.str. Germany 

Regensburg 

Donau 

Station 1 in Scheepmaker & Van Dalfsen (1989) 


g 

M 

G. fossarum s.str. Germany 

Regensburg 

Donau 

Station 2 in Scheepmaker & Van Dalfsen (1989) 


h 

M 

G. fossarum s.str. Germany 

Regensburg 

Donau 

Station 3 in Scheepmaker & Van Dalfsen (1989) 


i 

L 

G. fossarum s.l. 

Switzerland 

Oberaargau 

Langete 

Station 4 in Scheepmaker & Van Dalfsen (1989) 


j 

L 

G. fossarum s.l. 
G. p.pulex 

Switzerland 

Bern 

Aare 

Station 5 in Scheepmaker & Van Dalfsen (1989) 


k 

L 

G. fossarum s.l. 

Switzerland 

Bern 

Aare 

Station 6 in Scheepmaker & Van Dalfsen (1989) 


1 

K 

G. fossarum s.l. 

Belgium 

Luxembourg 

Meuse 

Station 7 in Scheepmaker & Van Dalfsen (1989) 


m 

K 

G. fossarum s.l. 

Belgium 

Lifcge 

Meuse 

Station 8 in Scheepmaker & Van Dalfsen (1989) 


n 

- 

G. fossarum s.l. 
G. p. pulex 

Belgium 

Namur 

Escaut 

Station 9 in Scheepmaker & Van Dalfsen (1989) 


o 

- 

G. fossarum s.l 

Belgium 

Brabant 

IJzer 

Station 10 in Scheepmaker & Van Dalfsen (1989) 


P 

J 

G. fossarum s.l. 

France 

Pas-de-Calais Liane 

Station 11 in Scheepmaker & Van Dalfsen (1989) 


q 

J 

G. fossarum s.l. 

France 

Pas-de-Calais Wimereux 

Station 12 in Scheepmaker & Van Dalfsen (1989) 


r 

J 

G. fossarum s.l. 

France 

Pas-de-Calais Slack 

Station 13 in Scheepmaker & Van Dalfsen (1989) 


s 

J 

G. fossarum s.l. 

France 

Pas -de-Calias Slack 

Station 14 in Scheepmaker & Van Dalfsen (1989) 




G.p. pulex 





t 


G. p. pulex 

Netherlands 

Gelderland 

Hierdense 

beek 

Station 15 in Scheepmaker & Van Dalfsen (1989) 



forms were described as G. stupendus Pinkster, 
1983 and G. orinos Pinkster & Scholl, 1984. The 
taxonomic status of the third form (G. cf. fossarum 
in Goedmakers & Roux, 1975, and Scheepmaker et 
al., 1988) has not yet been resolved. All three forms 
are genetically distinct from one another and from 
G. ibericus (cf. Pinkster & Scholl, 1984; Scheep¬ 
maker et al., 1988). 

Genetic differentiation of G. fossarum and G. 
caparti Petre-Stroobants, 1980 with reference to G. 
pulex was studied by Scheepmaker & Van Dalfsen 
(1989). These authors proposed the term G. fossa¬ 
rum sensu lato for all populations hitherto identi¬ 
fied as G. fossarum (including G. caparti) and 
reserved G. fossarum sensu stricto for the genetical¬ 
ly distinct populations from the area surrounding 
the former type locality near Regensburg (Federal 
Republic of Germany). 

This study primarily aims at resolving levels of 
genetic divergence among populations of G. fossa¬ 


rum sensu lato (cf. Scheepmaker & Van Dalfsen, 
1989), G. orinos , G. stupendus , and G. cf. fossa¬ 
rum, with reference to the “classical” species G. 
pulex - including G. pulex gallicus (S. Karaman, 
1931) - and G. wautieri . In addition, genetic differ¬ 
entiation was studied in relation to indirectly esti¬ 
mated levels of gene flow. Secondly, the relation 
between morphological and genetic differentiation 
among populations of the above species was inves¬ 
tigated. 

Genetic variation at 20 presumptive gene loci was 
studied using starch gel enzyme electrophoresis. 
The evolutionary significance of the genetic varia¬ 
tion recorded could be evaluated in part with the 
results from hybridization experiments by Roux 
(1971), Goedmakers (1972), Goedmakers & Roux 
(1975), Pinkster (1983) and Scheepmaker (1987). In 
addition to these experiments, crosses have been 
carried out between four population samples of G. 
fossarum sensu lato (including G. caparti ), one of 






Fig. 1A-G. Study areas and sampling localities (cf. Table I). 


G. orinos and one of G. cf. fossarum . 

Material and methods 

Sampling and collection sites 

Sampling for electrophoretic studies and cross¬ 
breeding experiments has been carried out accord¬ 
ing to the procedures of Scheepmaker (1987). In ad¬ 
dition to sampling by Scheepmaker et al., 1988 
(Fig. 1 A, areas F, H and I; Fig. IF, station a; Table 
I, stations a-e) and Scheepmaker & Van Dalfsen, 
1989 (Fig. 1 A, areas J-M, stations i, j and t; Table 
I, stations f-t), G. fossarum s.l. y G. orinos y G. cf. 
fossarum , and G. p. pulex were sampled at 25 
further localities. Moreover, G. wautieri and G. p . 
gallicus were sampled at 5 and 3 localities, respec¬ 
tively. The collection sites, sampling dates and the 
species collected are listed in Table I and Fig. 1. For 
detailed maps of areas I, H and F-M in Fig. 1A is 


referred to Scheepmaker (1987), Scheepmaker et al. 
(1988), and Scheepmaker & Van Dalfsen (1989), 
respectively. 

For a discussion and figures of morphological 
variation of G. fossarum s.l. in areas J, K, L and 
M (Fig. 1 A) is referred to Scheepmaker & Van Dalf¬ 
sen (1989). G. fossarum s.l . from areas B and C 
(Fig. 1 A) differs from G. fossarum s.l. population 
samples from other areas primarily by the total ab¬ 
sence of calceoli on the A2. The setation on the 
peduncle segments and flagellum of the A2 is regu¬ 
lar, comb-like, often longer than the segments on 
which they are implanted (for figures see Karaman 
& Pinkster, 1977). Individuals from area F resem¬ 
ble those from area L (Fig. 1A); individuals from 
station 17 are rather similar to those of areas B and 
C, but most individuals bear calceoli; and individu¬ 
als from station 18 (Fig. 1, area I) resemble in¬ 
dividuals from area J (Fig. 1 A). The morphological 
variation of some characters (cf. Scheepmaker & 
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Table II. Morphological differentiation of selected variable characters in G. fossarum s.l. Legend: » setae much longer than the di¬ 
ameter of the segment on which they are implanted; >= setae longer or equal to the diameter of the segment on which they are implant¬ 
ed; <= > setae shorter, equal or longer than the diameter of the segment on which they are implanted; - setation poor; + setation 
moderate; + + setation dense. 


Station 

no. 

Length of the setae with reference 
to the diameter of the segments 
on which they are inplanted 



Density setation 



Setation 

P3 

curly 

Calceoli 

present 


A2 


A2 










peduncle 


flagel- 










segment 


lum 











5 


P3 

P4 

A2 

P3 

P4 

U3 



>= 

>= 

>= 

» 

>= 

+ 

++ 

. 

+ 

+ 

_ 


>= 

» 

>= 

» 

= 

+ 

+/++ 

- 

+ 

+ 

- 

■i 

>= 

» 

>= 

» 

SB 

+ 

+ /++ 

- 

+ 

+ 

- 

■1 

>= 

» 

>= 

» 

= 

++ 

++ 

+ 

-/+ 

+ 

- 

mmm 

>= 

» 

>= 

» 

= 

++ 

++ 

+ 

-/+ 

+ 

- 


>= 

» 

>= 

» 

>= 

++ 

++ 

+ 

-/+ 

+ 

- 


>= 

>= 

>= 

» 


+ 

++ 

-/+ 

-/+ 

+ 

-/+ 

eg | 

>= 

>= 

>= 

» 

= 

+ 

++ 

-/+ 

-/+ 

+ 

-/+ 


>= 

>= 

>= 

> 

<=> 

+ 

+ 

- 

- 

+ 

- 

u 

= 

s 

>= 

» 

= 

- 

+ 

+ 

-/+ 

- 

+ 


Van Dalfsen, 1989) of G. fossarum s.l. in popula¬ 
tion samples from areas B, L, F, and I (+ station 
17) is listed in Table II. 

Morphological differentiation of the related 
forms G. orinos , G. cf. fossarum and G. stupendus 
of areas F, H, and I was discussed by Pinkster* & 
Scholl (1984), Scheepmaker et al. (1988), Pinkster 
(1983) and Scheepmaker (1987). Populations of G. 
cf. fossarum from area G (Fig. 1 A) are morphologi¬ 
cally similar to those of area H; likewise, popula¬ 
tions of G. orinos from areas D and E are similar 
to those of area F. The most conspicuous difference 
between G. wautieri and G. fossarum s.l. is the 
length of the inner ramus of uropod 3 (75% of the 
length of the outer ramus for G. wautieri , vs. 50% 
for G. fossarum). In females of G. wautieri the 
setation of the peduncle segments 4 and 5 is shorter 
than in females of G. fossarum (shorter than the 
diameter of the segment on which they are implant¬ 
ed), and the setae are implanted in a more irregular 
way. However, the younger females in particular 


are difficult to distinguish from G. fossarum s.l. 

Morphological differentiation of G. p. pulex and 
G. p. gallicus vs. G. fossarum s.l. and G. wautieri 
was reviewed by Karaman & Pinkster (1977). The 
morphological differentiation of some taxonomi- 
cally important characters among these species and 
forms is summarized in Table III; for figures see 
Karaman & Pinkster (1977) and Scheepmaker et al. 
(1988). 

Stations 18-20 (Fig. 1A, area I) correspond to 
stations 13, 17 and 43 in Scheepmaker (1987), and 
therefore they are not figured on a detailed map. As 
this study involved only 4 enzyme loci, new popula¬ 
tion samples were taken and analyzed. Stations 
15-17 were sampled to evaluate inter-specific, 
intra-areal comparisons between on the one hand 
local G. fossarum s.l. and G. orinos from area F, 
and on the other G. fossarum s.l. and G. stupendus 
from area 1+ station 17. Stations 9-14 include the 
localities where Roux (1971) and Goedmakers 
(1972) collected individuals of G. fossarum and 
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G. wautieri for cross-breeding experiments (viz. my 
stations 9 and 13 respectively). 

Cross-breeding experiments 

Cross-breeding experiments followed Pinkster 
(1983) and were carried out with samples of (a) G. 
fossarum s.l. ( G. caparti) from the type locality of 
G. caparti (area K, station 1), G. fossarum s.l. from 
area J (station r), G. fossarum s.l. from area L (sta¬ 
tion k), G. orinos from the type locality (area F, sta¬ 
tion a), and G. cf. fossarum from area H (station 
b). All combinations of crosses were performed, ex¬ 
cept for crosses among samples of G. orinos and G. 
fossarum s.l. from stations a and r, which were car¬ 
ried out by Pinkster & Scholl (1984). 


Electrophoresis 

Electrophoresis and staining procedures were iden¬ 
tical to those of Scheepmaker et al. (1988). The fol¬ 
lowing enzymes (E.C. numbers between paren¬ 
theses) were assayed: ADA — Adenosine deaminase 
(3.5.4.4); ALP-Alkaline phosphatase (3.1.3.1); 
APK - Arginine phosphate kinase (2.7.3.3); EST - 
Esterase (3.1.1.1); GDH-Glutamate dehydro¬ 
genase (1.1.1.47); GOT-Glutamic oxaloacetic 
transaminase (2.6.1.1); GPI —Glucose phosphate 
isomerase (5.3.1.11); HK - Hexokinase (2.7.1.1); 
LAP-Leucine aminopeptidase (3.4.1.1); MDH- 
Malate dehydrogenase (1.1.1.37); ME-Malic en¬ 
zyme (1.1.1.40); MPI —Mannose phosphate iso¬ 
merase (5.3.1.8); PEP-Peptidase (3.4.11/13); 
6PGD-6 Phosphogluconate dehydrogenase 
(1.1.1.44); PK-Pyruvate kinase (2.7.1.40). 

Analysis of allozyme variation 

The genetic interpretation of the variation was in¬ 
ferential. Data were analyzed by the computer pro¬ 
gram BIOSYS-1 (Swofford & Selander, 1981). Ma¬ 
trices of Rogers’ (1972) genetic distance, and the 
genetic distance and -identity according to Nei 
(1972; D and / respectively) were calculated for all 


population samples, including the French samples 
studied by Scheepmaker et al., 1988 (stations a-e) 
and samples studied by Scheepmaker & Van Dalf- 
sen, 1989 (stations f-t). From these data, an UPG- 
MA dendrogram (Sneath & Sokal, 1973) and a dis¬ 
tance Wagner network (Farris, 1972) were con¬ 
structed (for sake of space, complete distance 
matrices are not presented here). 

To measure intra-specific differentiation, 0- 
values according to Weir & Cockerham (1984) were 
calculated with the aid of a Macintosh computer 
program produced by W.N. Ellis (1989). For com¬ 
parison, F sr values according to Wright (1931, 
1951,1965) were also calculated. For loci with mul¬ 
tiple alleles, a weighted average 0 (and F ST ) was 
calculated. Significance of the observed 0-values 
was tested with a chi-square heterogeneity test 
(Workman & Niswander, 1970, in Waples, 1987). 
Quantitative gene flow estimates b/mp ST and Nm $ 
were derived from the relationship F ST - 
1/(1 +4Nm) in which N is the effective population 
size and m the migration rate (Wright, 1931). As in¬ 
dices of gene flow based on /^-statistics may be in¬ 
fluenced by natural selection and therefore may 
constitute poor indicators of the genetic structure 
of populations, these indices were used in combina¬ 
tion with the methods of Slatkin (1981, 1985). 
These methods do not depend on the assumption of 
selective neutrality (Barton & Slatkin, 1987). Gene 
flow levels among populations and areas were also 
calculated from the average frequency of private al¬ 
leles, p (1), using Slatkin’s (1985) formula: 

- In ([/?(!)]+ 2.44) 

Nm = e °' 505 

This formula was derived empirically from simula¬ 
tion experiments where N = 25. To adjust for 
differences in sample size, a corrected Nm value, 
Nm *, is obtained by dividing the Nm estimate by 
the ratio of the actual average sample size to 25 
(Slatkin, 1985): 

Nm * = i n which AT t is the actual average 

N 1 

J> act 

sample size. 
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Table III. Morphological differentiation of some taxonomically important characters in G. orinos , G. cf. fossarum , G. stupendus , 
G. fossarum s.L, G. waufieri , G. p. pulex and G. p. gatlicus. 


species 

G. orinos 

G. cf. fossarum 

G. stupendus 

G. fossarum s.L 

G. wautieri 

G. p. pulex 

G. p. gallicus 

Setosity A2 

moderate 

moderate 

moderate 

little 

little 

dense 

little 

Flagellum swollen 

- 

- 

- 

- 

- 

+ 

+ 

Calceoli 

+ 

- 

. 

-/+ 

+ 

+ 

+ 

Setosity P5-7 

long 

moderate 

long /moderate 

little 

absent 

absent 

absent 

Setae longer than spines 

++ 

+ 

-/+/++ 

- 

- 

- 

- 

Uropod 3: endop./exop. 

1/2 

1/2 

1/2 

1/2 

3/4 

3/4 

3/4 

Urosomites elevated 







+ 


Conditional average frequency (CAF; Slatkin, 
1981) gene flow profiles were obtained by plotting 
the average values, p (/), of electromorphs having 
the same incidence (i/d, where i is the number of 
population samples an electromorph appears in, 
and d the total number of populations sampled), 
against their incidence or occupancy number (Slat- 
kin, 1981; Caccone, 1985; Caccone & Sbordoni, 
1987). These procedures have been applied to the 
complete dataset of a species and to different sub¬ 
sets of populations. CAF gene flow profiles and p 
(l)-values have not been calculated for subsets of 
less than 4 population samples. 


Results 

Twenty presumptive gene loci from 17 enzyme sys¬ 
tems were scorable in all populations. Electro¬ 
morph frequencies and relative mobilities of the en¬ 
zymes coded for by these loci are listed in Table IV. 
Electromorphs of the Ada , Alp-l y Alp-2 9 Apk , 
Gdh , 6Pgd and Lap loci were monomorphic and 
fixed for the same electromorph in all population 
samples. The number of isozyme loci and quater¬ 
nary structures resolved were consistent with the in¬ 
terpretation of Scheepmaker et al. (1988). 

Genetic variation among populations 

Differences in electromorph distribution among the 
populations sampled for this study are shown in 


Table IV. The Est-l c and Got-1* electromorph 
almost completely discriminated G. p . pulex and G. 
p . gallicus from all other samples involved. G. p. 
pulex and G. p . gallicus were also characterized by 
the occurrence of the Mpi d and in most cases the 
Pk b electromorphs, and the Pep-2 C electromorph 
was only found in G. p . pulex. 

G. wautieri is characterized in most comparisons 
by the Hkf (also characterizing G. fossarum s.L 
from area F), Mpi 1 (shared with G. fossarum s.L 
from area F and station 17), PkA, Pep-2 m (except 
for G. wautieri from station 14), and Pep-4 h elec¬ 
tromorphs. G. orinos from area F is characterized 
by the monomorphic Mpi h electromorph. The 
Est-l e electromorph completely discriminates G. 
cf. fossarum from area G from all other samples 
(including G. ci. fossarum from area H); these sam¬ 
ples (and G. cf. fossarum from area H; Scheepma¬ 
ker et al., 1988) are also characterized by the Gpi k 
electromorph. The occurrence of the Gpi* electro¬ 
morph is limited to G. wautieri (areas B, C) and one 
sample of G. stupendus (station 19, area I). Local 
populations may deviate from the general pattern 
(thus contributing to a discontinuous local varia¬ 
tion pattern) by the occurrence of typical and highly 
frequent electromorphs (e.g. Pep-2 l and Mdh-l e in 
G. fossarum s.L from station 18, area I; Table I, 
Fig. 1 A; Got-1 8 and Gpi 8 in G. stupendus from sta¬ 
tion 19, area I). 

Differences between numbers of loci that are dis¬ 
criminative or diagnostic (0.99 criterion; Ayala & 
Powell, 1972) among pairwise comparisons of spe¬ 
cies or forms by area (including data from Scheep- 
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M. Scheepmaker - Gene flow in members of the Gammarus pulex-group 


Table IV. Electromorph frequencies at 13 presumptive gene loci, o = G. orinos ; cf. / = G. cf .fossarum; f = G. fossarum; s = G. 
stupendus; w = G. waufieri; g = G. p. gallicus; p = G. p. pulex ; N = sample size; h = heterozygosity per locus (direct count); H 
= mean heterozygosity over all (including monomorphic) loci; * = significant departure from Hardy-Weinberg distribution (P < 0.05); 
electromorph designation from fastest to slowest in alphabetical order, in accordance with Scheepmaker & Van Dalfsen (1989); in paren¬ 
theses: alternative allelic designation employed by Scheepmaker et al., 1988. 


station 


1 

2 

3 

4 

5 

6 

1 

8 

9 

10 

11 

12 

13 

14 

15 

16 



o 

o 

o 

0 

0 

o 

cf./ 

cf. / 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

locus 


















electromorph 
















Est-1 

(N) 

25 

24 

33 

30 

20 

25 

14 

20 

20 

23 

12 

15 

29 

22 

23 

26 


b 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

(b) 

c 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

(c) 

d 

1.00 

1.00 

1.00 

1.00 

0.75 

1.00 

0.00 

0.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 


e 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


f 

0.00 

0.00 

0.00 

0.00 

0.25 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


h 

- 

- 

- 

- 

0.30 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

Est-2 

(N) 

25 

24 

33 

30 

20 

25 

14 

20 

20 

23 

12 

15 

29 

22 

23 

26 


a 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

(a) 

b 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.28 

0.00 

0.00 

1.00 

0.00 

0.76 

0.92 

(b) 

c 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.72 

1.00 

1.00 

0.00 

1.00 

0.24 

0.08 


h 

- 

- 

- 

- 

- 

- 

- 

- 

- 

0.39 

- 

- 

- 

- 

0.30 

0.08 

Got-1 

(N) 

19 

24 

49 

26 

30 

33 

27 

23 

46 

45 

23 

31 

32 

30 

46 

33 


a 

0.00 

0.00 

1.00 

1.00 

0.95 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.90 

1.00 

(a) 

b 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.76 

1.00 

1.00 

0.98 

0.95 

1.00 

0.10 

0.00 


d 

0.00 

0.46 

0.00 

0.00 

0.05 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

(c) 

e 

1.00 

0.54 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


f 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


g 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.24 

0.00 

0.00 

0.02 

0.05 

0.00 

0.00 

0.00 


h 

- 

0.50 

- 

- 

0.03 

- 

- 

- 

0.43 

- 

- 

0.03 

0.09 

- 

0.02 

- 

Got-2 

(N) 

19 

24 

49 

26 

30 

33 

27 

23 

36 

48 

23 

31 

32 

30 

46 

32 

(a) 

b 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.99 

0.98 

1.00 

0.00 

0.00 

0.00 

0.00 

(b) 

c 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.00 

0.02 

0.00 

1.00 

1.00 

1.00 

0.98 


e 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.02 


h 

- 

- 

- 

- 

- 

- 

- 

- 

- 

0.02 

0.04 

- 

- 

- 

- 

0.03 

Gpi 

(N) 

32 

31 

49 

39 

41 

43 

37 

49 

62 

47 

28 

62 

42* 

47 

54 

26 


a 

0.20 

0.00 

0.00 

0.00 

0.04 

0.06 

0.00 

0.00 

0.77 

0.20 

0.77 

0.43 

0.24 

0.09 

0.00 

0.00 


b 

0.75 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.22 

0.11 

0.21 

0.01 

0.07 

0.80 

0.02 

0.00 


c 

0.00 

0.00 

0.99 

1.00 

0.61 

0.51 

0.00 

0.00 

0.01 

0.69 

0.02 

0.56 

0.69 

0.11 

0.92 

0.96 


d 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.01 

0.00 

(d) 

e 

0.05 

0.00 

0.01 

0.00 

0.29 

0.41 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.05 

0.02 

(e) 

g 

0.00 

0.00 

0.00 

0.00 

0.06 

0.02 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.02 


h 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


i 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


j 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


k 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


1 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


r 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


h 

0.31 

- 

0.02 

- 

0. 49 

0.53 

- 

- 

0.29 

0.34 

0. 39 

0.32 

0. 48 

0.25 

0.13 

0.08 

Hk 

(N) 

50 

38 

23 

40 

40 

43 

28 

33 

52 

35 

16 

29 

28 

11 

43 

18 

(a) 

b 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

(b) 

c 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

Mdh-1 

(N) 

20 

24 

54 

26 

30 

39 

27 

23 

46 

52 

26 

31 

32 

30 

39 

32 

(a) 

b 

0.45 

0.56 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.98 

0.97 

1.00 

1.00 

1.00 

1.00 

0.95 

1.00 


d 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

(b) 

e 

0.55 

0.44 

0.99 

0.02 

1.00 

0.83 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.05 

0.00 

(c) 

f 

0.00 

0.00 

0.01 

0.98 

0.00 

0.17 

0.00 

0.00 

0.02 

0.03 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


g 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


h 

0.40 

0. 71 

0.02 

0.04 

- 

0.23 

- 

- 

0.04 

0.06 

- 

- 

- 

- 

0.02 

- 




Bijdragen tot de Dierkunde, 60 (1) - 1990 
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Table IV. Cont. 


station 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 



o 

o 

o 

0 

0 

o 

cf.f 

cf./ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

locus 


















electromorph 
















Mdh-2 

(N) 

19 

24 

54 

26 

30 

33 

27 

23 

36 

45 

23 

31 

32 

30 

45 

31 


a 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.03 

0.00 

0.00 

0.00 

0.00 

0.00 

0.01 

0.00 


b 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


c 

0.97 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.97 

1.00 

1.00 

1.00 

1.00 

1.00 

0.99 

1.00 


e 

0.03 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


h 

0.05 

- 

“ 

- 

- 

- 

- 

- 

0. 07 

- 

- 

- 

- 

- 

0. 20 

- 

Me 

(N) 

28 

37 

47 

61 

43 

41 

36 

39 

50 

56 

25 

14 

43 

35 

33 

30 


b 

1.00 

1.00 

1.00 

0.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


c 

0.00 

0.00 

0.00 

1.00 

0.00 

0.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 


d 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


e 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

Mpi 

(N) 

27 

31 

52 

39 

44 

35 

37 

40 

62 

55 

23 

45* 

41 

46 

62 

26 


b 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


d 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


e 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

(b) 

f 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

h 

1.00 

1.00 

0.02 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


i 

0.00 

0.00 

0.00 

0.00 

0.49 

0.39 

0.00 

0.00 

0.00 

0.06 

0.00 

0.00 

0.00 

0.00 

0.97 

1.00 


j 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.02 

0.00 

0.00 

(0 

k 

0.00 

0.00 

0.52 

1.00 

0.00 

0.00 

0.00 

0.04 

0.04 

0.05 

0.00 

0.04 

0.05 

0.03 

0.03 

0.00 

(g> 

m 

0.00 

0.00 

0.46 

0.00 

0.51 

0.60 

1.00 

0.95 

0.94 

0.30 

0.24 

0.34 

0.93 

0.94 

0.00 

0.00 

n 

0.00 

0.00 

0.00 

0.00 

0.00 

0.01 

0.00 

0.01 

0.00 

0.59 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.02 

0.00 

0.76 

0.62 

0.02 

0.01 

0.00 

0.00 


h 

- 

- 

0. 44 

- 

0. 52 

0. 51 

- 

0.1 

0.13 

0.51 

0.22 

0.27 

0.15 

0. 13 

0.06 

- 

Pep-2 

(N) 

29 

27 

29* 

39 

40* 

40* 

28 

31 

32 

39 

23 

24 

16 

28 

26 

41 


c 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


d 

0.00 

0.00 

0.00 

0.00 

0.18 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


e 

0.00 

0.00 

0.78 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

(e) 

f 

0.00 

0.00 

0.00 

0.00 

0.12 

0.26 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

h 

1.00 

0.00 

0.00 

0.00 

0.70 

0.74 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


j 

0.00 

1.00 

0.22 

0.00 

0.00 

0.00 

0.00 

0.00 

0.80 

1.00 

0.28 

0.94 

1.00 

0.98 

1.00 

1.00 


k 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.20 

0.00 

0.72 

0.06 

0.00 

0.02 

0.00 

0.00 


1 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


m 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


n 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


h 

- 

- 

0.17 

- 

0. 25 

0.22 

- 

- 

0.22 

- 

0.22 

0.12 

- 

0. 04 

- 

- 

Pep-4 

(N) 

29 

27 

29 

39 

39 

41 

28 

31 

32 

39 

26 

25 

19 

28 

26 

41 

(a) 

c 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.04 

0.00 

d 

1.00 

1.00 

1.00 

1.00 

0.89 

0.92 

1.00 

0.95 

1.00 

0.90 

1.00 

0.58 

0.92 

1.00 

0.96 

0.92 


e 

0.00 

0.00 

0.00 

0.00 

0.11 

0.08 

0.00 

0.05 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.08 


f 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


g 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.03 

0.00 

0.00 

0.00 


h 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.10 

0.00 

0.42 

0.05 

0.00 

0.00 

0.00 


h 

* 

- 

- 

- 

0.18 

0.12 

- 

0.10 

- 

0. 20 

- 

0.28 

0.10 

- 

0.08 

0.12 

Pk 

(N) 

40 

37 

49 

49 

33 

46 

30 

36 

36 

39 

23 

38 

20 

25 

26 

41 


b 

0.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

(b) 

c 

1.00 

0.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

d 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


e 

h 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


H 

0.04 

0.06 

0.03 

0.00 

0. 09 

0. 08 

0.00 

0. 01 

0.06 

0.08 

0.04 

0.06 

0. 04 

0. 02 

0.04 

0. 01 
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M. Scheepmaker - Gene flow in members of the Gammarus pulex-group 


Table IV. Cont. 


station 17 18 19 20 11 12 13 14 21 22 23 24 23 

ffsswwwwwgggp 

locus _ 

electromorph 


Est-1 


(N) 

29 

31 

30 

24 

11 

5 

31 

20 

18 

28 

35 

25 

10 

b 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

c 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 

1.00 

d 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

e 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

f 

h 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


Est-2 

(a) 

(b) 

Got-1 

(a) 

(c) 


Got-2 

(a) 

(b) 

Gpi 


(d) 

(e) 


Hk 

(a) 

(b) 

Mdh-1 

(a) 

(b) 

(c) 


(N) 

29 

31 

30 

24 

11 

5 

31 

20 

18 

28 

35 

25 

10 

a 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

b 

0.00 

1.00 

0.02 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

c 

1.00 

0.00 

0.98 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

h 

- 

- 

0.03 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 


(N) 

18 

12 

20 

7 

3 

9 

27 

9 

26 

27 

25 

28 

18 

a 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

0.98 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

b 

1.00 

1.00 

0.00 

0.71 

0.00 

0.00 

0.02 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

d 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

c 

0.00 

0.00 

0.00 

0.29 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

f 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

0.88 

1.00 

1.00 

g 

0.00 

0.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.12 

0.00 

0.00 

h 

- 

- 

- 

0. 29 

- 

- 

0.04 

- 

- 

- 

0.24 

- 

- 


(N) 

18 

12 

20 

9 

3 

9 

27 

9 

26 

28 

25 

28 

18 

b 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

c 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

e 

h 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

(N) 

33 

30 

73 

39 

20 

15 

31 

30 

59 

43 

44 

34 

18 

a 

0.89 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

b 

0.02 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.01 

0.00 

0.00 

c 

0.09 

0.00 

0.00 

0.47 

0.00 

0.00 

0.00 

0.00 

0.00 

0.32 

0.42 

0.82 

0.89 

d 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

e 

0.00 

0.00 

0.02 

0.44 

0.00 

0.00 

0.00 

0.00 

0.00 

0.04 

0.36 

0.07 

0.00 

g 

0.00 

0.00 

0.96 

0.06 

0.65 

0.00 

1.00 

0.95 

0.97 

0.64 

0.21 

0.00 

0.06 

h 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.10 

0.00 

i 

0.00 

0.00 

0.00 

0.03 

0.35 

1.00 

0.00 

0.05 

0.03 

0.00 

0.00 

0.00 

0.00 

j 

0.00 

0.00 

0.01 

0.00 

0.00 

0.00 

0,00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

k 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1 

0.00 

0.00 

0.01 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

r 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.05 

h 

0. IS 

- 

0.10 

0.38 

0. 20 

- 

- 

0.10 

0. 07 

0.42 

0. 57 

0. 24 

0. 22 


(N) 

39 

21 

34 

39 

15 

2 

25 

18 

41 

33 

20 

23 

13 

b 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 

1.00 

c 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 


(N) 

18 

11 

20 

9 

3 

9 

28 

9 

26 

28 

25 

38 

16 

b 

1.00 

0.00 

1.00 

1.00 

0.00 

1.00 

0.91 

1.00 

1.00 

0.93 

0.84 

0.97 

0,53 

d 

0.00 

0.00 

0.00 

0.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.07 

0.00 

0.00 

0.00 

e 

0.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.09 

0.00 

0.00 

0.00 

0.16 

0.03 

0.47 

f 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

g 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

h 

- 

- 

- 

- 

- 

- 

0. 18 

- 

- 

0.14 

0. 24 

0.05 

0. 31 
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Table IV. cont. 


station 


17 

18 

19 

20 

11 

12 

13 

14 

21 

22 

23 

24 

25 



/ 

/ 

s 

s 

w 

w 

w 

w 

w 

8 

8 

8 

P 

locus 















electromorph 













Mdh-2 

(N) 

18 

11 

20 

9 

3 

9 

28 

9 

26 

28 

25 

38 

18 


a 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


b 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


c 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 


e 

h 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

Me 

(N) 

26 

46 

54 

47 

7 

6 

37 

19 

31 

25 

30 

20 

11 


b 

0.00 

0.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 


c 

0.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 

1.00 


d 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


e 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

Mpi 

(N) 

30 

47 

43 

40* 

13 

10 

24 

19 

47 

48* 

18 

23 

13 


b 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.06 

0.96 

0.00 


d 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.84 

0.80 

0.04 

0.92 


e 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.01 

0.11 

0.00 

0.00 

(b) 

f 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.15 

0.00 

0.00 

0.08 

h 

0.00 

0.02 

0.45 

0.23 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.03 

0.00 

0.00 


i 

1.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 


j 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

(0 

k 

0.00 

0.00 

0.28 

0.41 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

(g) 

m 

0.00 

0.35 

0.27 

0.33 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

n 

0.00 

0.00 

0.00 

0.04 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


0 

0.00 

0.63 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


h 

- 

0.36 

0. 67 

0.35 

- 

- 

- 

- 

- 

0.08 

0.28 

0. 09 

0.15 

Pep-2 

(N) 

24 

23 

33 

29 

25 

6 

26 

28 

28 

31 

23 

18 

13 


c 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 


d 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

0.00 


e 

0.00 

0.00 

0.08 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

(e) 

f 

0.00 

0.00 

0.89 

0.85 

0.00 

0.00 

0.00 

0.00 

0.00 

0.55 

0.59 

0.00 

0.00 

h 

0.00 

0.00 

0.00 

0.15 

0.00 

0.00 

0.00 

0.00 

0.00 

0.31 

0.28 

0.00 

0.00 


j 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.07 

0.14 

0.13 

0.00 

0.00 


k 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.52 

0.00 

0.00 

0.00 

0.00 

0.00 


1 

0.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


m 

0.00 

0.00 

0.03 

0.00 

1.00 

1.00 

1.00 

0.00 

0.93 

0.00 

0.00 

0.00 

0.00 


n 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.48 

0.00 

0.00 

0.00 

0.00 

0.00 


h 

’ 

- 

0. 21 

0.24 

- 

- 

- 

0.32 

0. 07 

0.35 

0.48 

- 

- 

Pep-4 

(N) 

16 

23 

33 

30 

25 

6 

24 

28 

28 

31 

23 

18* 

13 

(a) 

c 

0.00 

0.00 

0.94 

0.95 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.02 

0.17 

0.00 

d 

1.00 

1.00 

0.06 

0.05 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

0.85 

0.83 

0.92 


e 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


f 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.08 


g 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


h 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.00 

0.13 

0.00 

0.00 


h 

- 

- 

0.12 

0.10 

- 

- 

- 

- 

- 

- 

0. 30 

0.00 

0.15 

Pk 

(N) 

14 

41 

45 

43 

25 

2 

25 

28 

29 

51 

28 

38 

23 


b 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.97 

1.00 

1.00 

1.00 

(b) 

c 

1.00 

1.00 

0.94 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.02 

0.00 

0.00 

0.00 

d 

0.00 

0.00 

0.06 

0.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 


e 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


h 

- 

- 

0. 07 

- 

- 

- 

- 

- 

- 

0.06 

- 

- 

- 


H 0.01 0.02 0.06 0.07 0.01 


0.01 0.02 0.01 0.05 0.10 0.02 0.04 
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M. Scheepmaker - Gene flow in members of the Gammarus pulex-group 


Table V. Number of diagnostic loci (0.99 criterion; Ayala & Powell, 1972) among species and forms, averaged by area. Data from 
Scheepmaker et al. (1988) and Scheepmaker & Van Dalfsen (1989) are included. 


area(sy 

station(s) no. of 
samples 


L B+C F 1+17 K F DfE H 


I B+C J+L I+H 
25, j, t 


G. fossarum 

J 

5 

_ 











s. 1. 














G. fossarum 

j 

l 

2 

- 










s. 1. 














G. fossarum 

M 

3 

4 

2* 

- 









s. str. 














G. fossarum 

L 

3 

1 

- 

2 

- 








s. /. 














G. fossarum 

B+C 

6 

- 

2 

4 

- 

- 







s. /. 














G. fossarum 

F 

2 

1 

3 

4 

1 

1 

- 






s. 1. 














G. fossarum 

1+17 

2 

1 

3 

4 

- 

- 

1 

- 





s. 1. 














G. fossarum 

K 

2 

- 

2 

2 

- 

- 

1 

- 

- 




s. L 














G. orinos 

F 

3 

3 

2 

4 

2 

4 

4 

2 

2 

- 



G. orinos 

DfE 

4 

- 

- 

2 

- 

- 

- 

- 

- 

- 

- 


G. cf. fossarum 

H 

3 

5 

4 

5 

3 

5 

4 

4 

5 

4 

1 

- 

G. cf. fossarum 

G 

2 

4 

5 

6 

3 

5 

5 

4 

4 

4 

2 

4 

G. stupendus 

I 

3 

1 

- 

2 

- 

2 

3 

2 

2 

- 

- 

3 3 

G. wautieri 

B+C 

5 

6 

7 

4 

3 

6 

3 

6 

6 

6 

4 

5 7 5 - 

G. p. pulex 

J+L, 

5 

3 

4 

5 

4 

5 

6 

5 

4 

4 

3 

7 4 4 6 - 


25, j, t 
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maker et al., 

1988 and Scheepmaker & Van Dalf- 


0.63 

-0.69 t D 


0.47-0.38). As a consequence, 


sen, 1989; Table I, stations a-t) are summarized in 
Table V. Values of /and D , averaged by taxon, and 
by area, are listed in Tables VI and VII, respective¬ 
ly. Because of the low levels of divergence recorded, 
and the considerable geographic distances separat¬ 
ing most sampling localities (Table I; Fig. 1A), G . 
p . pulex and G. p. gallicus were averaged by form, 
not by area. In the UPGMA dendrogram (Fig. 2) 
and the Wagner distance network (Fig. 3), clearly 
distinct clusters are formed by G. fossarum s.str., 
G. stupendus, G. wautieri and G. p . pulex + G. p . 
gallicus. Intra-specific inter-areal divergence may 
be considerable: inter-areal /-values among G. ori- 
nos from areas D/E, and F range from 0.67-0.87 
(D = 0.40-0.14); among G. cf. fossarum from 
areas G and H inter-areal /-values range from 


neither the conspecific G. orinos samples from 
areas D/E, and F, nor the conspecific G. ci. fossa¬ 
rum samples from areas G and H form mono- 
phyletic groups in Figs. 2 and 3. 

Inter-specific inter-areal /-values among G. ori¬ 
nos and G. cf. fossarum do not differ significantly 
from intra-specific inter-areal values in a Mann- 
Whittney U-test at the 0.05 level (except for G. ori¬ 
nos from areas D, E and F, vs. G. cf. fossarum 
from area H and G. orinos from area F). Moreover, 
inter-areal /-values between any of the combina¬ 
tions F, D/E, G and H do not differ significantly 
at the 0.05 level. To put it another way, G. orinos 
from areas D/E seems to differ about as much from 
G. orinos from area F, as it does from G. cf./ass<ar- 
rum from area G, or H, and vice versa. The same 


Table VI. Matrix of Nei’s (1972) genetic identity coefficients among species and forms averaged by area, based upon 20 presumptive enzyme loci. Data from Scheepmaker 
et al. (1988) and Scheepmaker & Van Dalfsen (1989) are included. Range / m j n -/ max in italics. 
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G.wautieri B+C 5 0.60 0.66 0.68 0.63 0.58 0.71 0.58 0.60 0.64 0.67 0.64 0.51 0.68 0.95 

(0.55-0.62) (0.63-0.68) (0.64-0.71) (059-0.65) (052-0.63) (0.67-0.73) (051-0.66) 056-0.62) (0.62-0.67) (0.60-0.72) (0.63-0.66) (050-051) (0.62-0.73) (0.91-1.00) 
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Fig. 2. UPGMA dendrogram of Nei’s (1972) genetic distance ( D ) based upon 20 enzyme loci. 
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Fig. 3. Distance Wagner network of Rogers’ (1972) genetic distance rooted at the midpoint of the longest path (Farris, 1972) based upon 
20 enzyme loci. 





Table VII. Matrix of Nei’s (1972) genetic distance coefficients among species and forms averaged by area, based upon 20 presumptive enzyme loci. Data from Scheepmaker 
et al. (1988) and Scheepmaker & Van Datfsen (1989) are included. Range D min -D max in italics. 
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Table VIII. 0-indices according to Weir & Cockerham (1984) of the species examined. Data from Scheepmaker et al. (1988), and Scheep- 
maker & Van Dalfsen (1989) are included. G. foss. s.L = G. fossarum s.i G. foss. s. sir. = G. fossarum s. sir.; G. or. = G. orinos; 
G. cf .foss. = G. cf. fossarum; G. stup. - G. stupendus; G. waul. = G. wautieri ; G. p.p. = G. p. pw/e*; G. p. g. = G. p. gallicus. 
Significance of population differentiation was tested with a chi-square heterogeneity test; * P < 0.05. 


locus 

G. foss. s. 1. 

G. foss. s. str. 

G. or. 

G. cf .foss. 

G. stup. 

G. waut. 

G. p. p. 

G. p. g. 

Apk 

_ 

0.15* 

_ 

_ 

_ 


_ 

_ 

Est-I 

- 

0.36* 

- 

1.00* 

- 

- 

- 

- 

Est-2 

0.63* 

- 

- 

- 

0.01 

- 

- 

- 

Got-I 

0.56* 

0.04* 

0.90* 

1.00* 

0.89* 

0.02 

- 

0.11* 

Got-2 

0.77* 

0.25* 

- 

- 

- 

- 

0.00 

- 

Gpi 

0.57* 

0.20* 

0.58* 

0.50* 

0.83* 

0.22* 

0.06* 

0.28* 

Hk 

1.00* 

- 

- 

- 

- 

- 

- 

- 

Mdh-1 

0.30* 

0.06* 

0.72* 

- 

0.07* 

0.63* 

0.39* 

0.07* 

Mdh-2 

0.00 

0.01 

0.01 

0.40* 

- 

- 

- 

- 

Me 

1.00* 

- 

1.00* 

- 

- 

- 

- 

- 

Mpi 

0.59* 

0.46* 

0.72* 

0.96* 

0.40* 

- 

0.49* 

0.67* 

Pep-2 

0.73* 

0.05* 

0.76* 

1.00* 

0.15* 

0.68* 

0.02 

0.47* 

Pep-4 

0.34* 

- 

0.05* 

0.97* 

0.05* 

- 

0.38* 

0.10* 

Pk 

0.13* 

- 

1.00* 

- 

0.06 

- 

- 

- 

Mean 

0.63 

0.20 

0.75 

0.86 

0.67 

0.59 

0.38 

0.42 


Table IX. 0-indices according to Weir & Cockerham (1984) of subsets of the species examined. For legend see Table VIII. 


Area(s) 

Stations 

G. foss. s.L G. foss. s.L 
B, C, J, L B, C 

9-14; k-m; 9-14 

o-s 

G. foss. sJ. G. foss. s.l. G. foss. s.L G. foss. s.L 

C J J K, L 

9-12; 14 o-s p-s i;k-m 

G. or. 

D, E 

3-6 

G. or. 

F 

a; 1-2 

G. cf. foss. 

H 

b-d 

G. waut. 

B, C 
12-13; 15 

locus 

Est-2 

0.87* 

0.59* 

0.53* 




0.25* 




Got-I 

0.25* 

0.17* 

0.20* 

0.00 

0.01 

0.22* 

0.03 

0.52* 

- 

0.01 

Got-2 

0.76* 

0.98* 

0.98* 

0.04 

0.04 

0.45* 

- 

- 

- 

- 

Gpi 

0.53* 

0.40* 

0.42* 

0.35* 

0.00 

0.16* 

0.27* 

0.12* 

0.49* 

0.01 

Mdh-1 

0.19* 

0.00 

0.00 

0.00 

0.00 

0.14* 

0.80* 

0.28* 

- 

0.06 

Mdh-2 

0.01 

0.01 

0.01 

- 

- 

0.01 

- 

- 

0.06* 

- 

Mpi 

0.59* 

0.49* 

0.50* 

0.15* 

0.19* 

0.10* 

0.47* 

- 

0.00 

- 

Pep-2 

0.50* 

0.48* 

0.48* 

0.02 

0.05 

- 

0.60* 

1.00* 

0.02 

0.03 

Pep-4 

0.36* 

0.26* 

0.29* 

0.00 

0.00 

0.25 

0.05** 

0.01 

0.38* 

- 

Pk 

0.13* 

- 

- 

- 

- 

• 

• 

- 

■ 

* 

Mean 

0.48 

0.59 

0.53 

0.18 

0.12 

0.25 

0.60 

0.56 

0.31 

0.03 


holds true for each of these subsets with regard to 
G. stupendus and G. fossarum s.L . Although inter¬ 
areal /-values of this combination are not statisti¬ 
cally significant in a Mann-Whittney U-test, the 
highest intra-specific inter-areal value of /(/= 0.87) 


was found among G. orinos from areas D, E and F. 

G. fossarum s.L from areas I, F, B/C, L, K and 
J, plus stations 17 and o (Fig. 1A, Table I), 
although genetically well differentiated (/ m j n -/ max 
= 0.69-0.99), forms a delimited monophyletic 
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Table X. 0-indices according to Weir & Cockerham (1984) of pairwise combinations of samples of the species examined. For legend 
see Table VIII. 


Area(s) 

Stations 

G. foss. s.l. G. foss. s.l. 
C C 

10, 12 9, 14 

G. foss. sJ. G. foss. s.l. 
K L 

1, m i, k 

G. foss. s.1. 
F 

15, 16 

G. foss. s.l. 

I 

17, 18 

G. or. 

D, E 

5.6 

G. or. 

F 

a; 1 

G. cf. foss. 
H 
c, d 

G. cf. foss. 
G 

7,8 

G. p. g. 

G 

22, 23 

locus 

Est-1 







0.25* 





Est-2 

0.22* 

- 

- 

. 

0.07 

1.00* 

- 

- 

- 

- 

_ 

Got-l 

0.01 

0.20* 

0.05* 

0.13* 

0.07* 

- 

0.03 

0.01 

- 


0.11 

Got-2 

0.00 

- 

0.02* 

0.38* 

0.01 

- 

- 

- 

- 

- 

. 

Gpi 

0.06* 

0.54* 

0.00 

0.05* 

0.00 

0.90* 

0.01 

0.00 

0.01 

- 

0.20 

Mdh-1 

0.02 

0.00 

0.12* 

0.08* 

0.03 

1.00* 

0.07* 

0.42* 

- 

- 

0.07 

Mdh-2 

- 

0.01 

0.00 

- 

0.00 

- 

- 

0.03 

0.04 

. 

- 

Mpi 

0.41* 

0.41* 

0.06* 

0.12* 

0.01 

0.72* 

- 

- 

- 

0.03 

0.03 

Pep-2 

0.07 

0.13* 

- 

- 

- 

1.00* 

0.04 

1.00* 

- 


0.02 

Pep-4 

0.23* 

- 

0.49* 

0.05* 

0.02 

- 

0.01 

0.01 

- 

0.03 

0.14 

Pk 

* 


- 

- 

- 

- 

- 

■ 

- 

- 

0.14 

Mean 

0.25 

0.35 

0.29 

0.15 

0.04 

0.93 

0.05 

0.61 

0.02 

0.03 

0.09 


group of genetically related forms in Figs. 2 and 3. 
The recorded variation is not always congruent 
with the geographic distribution and appears to be 
gradual; sharp discontinuities are often lacking. 
However, due to increased local variation, popula¬ 
tion samples from areas B/C in Figs. 2 and 3 do not 
cluster together and they are occasionally closer to 
population samples from other areas (e.g. G. fossa- 
rum s.l. from stations 9 and 14 with G. fossarum 
s.l. from area J). G. fossarum s.l. from stations 17 
and 18 (area I) are genetically rather distant and 
cluster apart in Figs. 2 and 3. However, G. fossa¬ 
rum s.l. from stations j and n deviate strongly from 
this pattern; in Fig. 2 they cluster close to G. 
stupendus and G. orinos , and in Fig. 3 to G. pulex 
and to G. fossarum s. str. f respectively. 

Genetic structure of populations and levels of gene 
flow 

The results of the /^-statistic analysis (0) for G. fos¬ 
sarum s.l ., G. p. pulex , G. p. gallicus , G. wautieri , 
G. orinos , G. stupendus and G. cf. fossarum are 
shown in Table VIII. In the following procedures, 
because of their strongly deviating character, G. 
fossarum s.l. from stations j and n have been omit¬ 
ted. Significant population subdivision was found 


for all species at a high number of polymorphic loci 
(e.g. all 7 in G. cf. fossarum ; 8 of 9 in G. orinos ; 
11 of 12 in G. fossarum s.l.). Creating subsets by 
omitting genetically deviating samples or by com¬ 
bining stations by area (Table IX) does not change 
this pattern in most cases; however, much less 
population subdivision was found in subsets 
formed by G. fossarum s.l. (2 of 7 loci in stations 
o-s, and 1 of 7 in stations p-s) and G. wautieri 
(none in stations 12, 13 and 15). Subsets formed by 
pairwise combinations of genetically similar and 
geographically close populations (Table X) show 
that within several species and areas statistically sig¬ 
nificant population subdivision still exists over 
short distances (e.g. stations i-k, all 6 loci; stations 
1-m, 5 of 7 loci; stations 10-12 and 9-14, 4 of 8, 
and 4 of 6 loci, respectively). 

Mean values of 0 and F ST within species and sub¬ 
sets of population samples are listed in Tables 
XI-XIII, along with the corresponding Nm- 
values. For sets of 4 population samples and more, 
the average frequencies of private alleles p (1) and 
the corresponding Nm- and Nm*-v alues figure in 
Tables XI-XII. Values of 0 are generally higher 
than those of F sr Consequently, levels of gene 
flow derived from 0 are lower than those derived 
from f st . Genetic drift and bottlenecks will result 
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Table XI. Various Nm estimates in G. fossarum s.l. , G. fossarum s. str. t G. wautieri , G. stupendus , G. or in os, G. cf. fossa rum, G. 
p. pulex and G. p. gallicus. p (I) = average frequency of private alleles (Slatkin, 1985); Nav = average sample size; Nm = estimate 
derived from p (1); Nm* = idem, with correction for sample size; F $T - F-statistic (Wright, 1931); 0 = idem (Weir & Cockerham, 
1984); Nm FsT = Nm estimate derived from F ST \ Nm e = idem, derived from 0; n.p.a. = number of private alleles; km (av) = average 
geographic distance among samples in kilometres. 


Species 

Area(s) 

Stations 

No. of 
samples 

km (av) 

F sr 

9 

n.p.a 

m 

Nav 

Nm 

Nm* 

NmFsT 

Nm g 

G. fossarum sJ. 

B, C, F, I-L 

9-18; j-m; o-s 

19 

306 

0.68 

0.63 

12 

0.19 

33.31 

0.21 

0.16 

0.12 

0.15 

G. fossarum s. str. 

M 

f-h 

3 

25 

0.16 

0.20 






1.31 

0.99 

G. wautieri 

B,C 

11-14; 21 

5 

87 

0.74 

0.59 

6 

0.36 

13.03 

0.06 

0.11 

0.09 

0.18 

G. stupendus 

I 

19,20; e 

3 

21 

0.47 

0.67 






0.28 

0.13 

G. orinos 

D, E, F 

1-6; a 

7 

93 

0.71 

0.75 

7 

0.35 

33.93 

0.06 

0.05 

0.10 

0.08 

G. cl. fossarum 

H, G 

7,8; b-d 

5 

39 

0.86 

0.86 

6 

0.18 

37.65 

0.25 

0.16 

0.04 

0.04 

G. p. pulex 

J, L 

25; j, n, s, t 

5 

352 

0.36 

0.38 

8 

0.11 

21.06 

0.64 

0.76 

0.44 

0.40 

G. p. gallicus 

G, I 

22-24 

3 

169 

0.35 

0.42 






0.46 

0.35 


Table XII . Various /Vm-estimates in subsets of G. fossarum s.l ., 

G. orinos , G. 

cf .fossarum, and G. p. pulex. For legend see Table XI. 

Species 

Area(s) 

Stations 

No. of 
samples 

km (av) 

F st 

9 

n.p.a 

m 

Nav 

Nm 

Nm* 

N m F sT 

Nm, 

G. fossarum s. /. 

B, C.J-L 

9-14; i, k-m; o-s 

15 

289 

0.50 

0.48 

8 

0.03 

34.89 

8.26 

5.92 

0.25 

0.27 

G. fossarum s. 1. 

B,C 

9-14 

6 

35 

0.55 

0.59 

5 

0.12 

29.94 

0.50 

0.42 

0.20 

0.17 

G. fossarum s. /. 

B, C 

9-12; 14 

5 

12 

0.50 

0.53 






0.26 

0.22 

G. fossarum s. 1. 

J 

o-s 

5 

37 

0.15 

0.18 

5 

0.06 

32.53 

1.96 

1.51 

1.41 

1.17 

G. fossarum s. 1. 

J 

p-s 

4 

8 

0.10 

0.12 






2.20 

1.75 

G. fossarum s. 1. 

K, L 

i, k-m 

4 

270 

0.21 

0.25 

9 

0.06 

45.58 

2.03 

1.11 

0.93 

0.75 

G. orinos 

D, E 

3-6 

4 

35 

0.53 

0.60 

6 

0.25 

34.10 

0.12 

0.09 

0.22 

0.17 

G. orinos 

F 

1,2; a 

3 

9 

0.55 

0,56 






0.20 

0.19 

G. cf. fossarum 

H 

b-d 

3 

13 

0.27 

0.31 






0.68 

0.56 

G. p. pulex 

J , L 

25; n, s, t 

4 

292 

0.32 

0.38 






0.53 

0.40 

G. wautieri 

B, C 

12,13; 15 • 

3 

35 

0.02 

0.03 






12.25 - 

8.08 


Table XIII. Nm estimates derived from F ST (Wright, 1931) and 
0 (Weir and Cockerham, 1984) in pairwise combinations of 
G. fossarum s.l. , G. orinos, G. cf. fossarum, and G. p. gallicus . 
For legend see Table VIII. 


Species 

Area(s) 

Stations 

km (av) F„ 

9 

Nm Fsr 

Nm e 

G, fossarum s. 1. 

C 

10, 12 

5 

0.15 

0.25 

1.39 

0.74 

G. fossarum s. 1. 

C 

9, 14 

4 

0.23 

0.35 

0.86 

0.47 

G. fossarum s. 1. 

I 

17. 18 

55 

0.89 

0.93 

0.03 

0.02 

G. fossarum s. L 

F 

15. 16 

3 

0.03 

0.04 

7.10 

5.56 

G. fossarum s. 1. 

K 

l.m 

40 

0.19 

0.29 

1.05 

0.63 

G. fossarum s. 1. 

L 

i,k 

41 

0.08 

0.15 

2.91 

1.42 

G. orinos 

D.E 

5, 6 

65 

0.03 

0.05 

7.10 

4.79 

G. orinos 

F 

1; a 

3 

0.42 

0.61 

0.35 

0.16 

G. cf. fossarum 

H 

c, d 

10 

0.02 

0.02 

15.38 - 

11.38- 

G. cf. fossarum 

G 

7, 8 

4 

0.02 

0.03 

10.62 - 

8.51 

G. p. gallicus 

G 

22, 23 

7 

0.06 

0.09 

4.21 

2.53 


in substantial local differentiation if Nm < l;when 
Nm > 1, populations tend to become panmictic for 
neutral alleles (Wright, 1931). In most species and 
subsets examined, Nm*, Nm FsT , and Nm e are < 1 
(Tables XI-XIII), and these values are much low¬ 
er than expected based on dispersal abilities in a 
stable aquatic environment. At first sight, the im¬ 
portant geographic distances separating some of 
the sampling localities of populations forming sets 
and subsets in Tables XI-XIII are likely to account 
for the low Nm-values observed. However, based 
upon geographic distance data in Tables XI-XIII, 
there does not seem to be a correlation with 
Nm e - values. 
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G. p. pulex 

(S-5) 


G. fossarum s.l. 

(S = 5; stations o-s) 

-► 

G. orinos 

(S-7) 

► 

G. fossarum s.l. 

(S - 6; stations 9-14) 

◄-► 

G. cf. fossarum 

(S-5) 

+ -► 

G. orinos 

(S - 4; stations 3-6) 


G. wautieri 

(S-5) 


G. fossarum s.l. 

(S - 15; stations 9-14, i, k-m, o-s) 

— 

G. fossarum s.l. 

(S - 19) 

— - 

G. fossarum s.l. 

(S = 6; stations i, K-m) 



Fig. 4. CAF gene flow profiles (Slatkin, 1981) in G. p. pulex , G. 
or in os, G. cf. fossa rum, G. wautieri and G. fossarum s.l. 
S = number of samples; p (/) = conditional average frequency 
(CAF) of an allele; i/d = incidence of an allele, where / = the 
total number of samples it appears in, and d = the total number 
of samples examined. The points plotted are average values of 
p (/) for all alleles having the same i/d value. 

Slatkin’s (1985) approach of private alleles 
(Nm*) has been applied to 10 sets of populations. 
Only once (Table XII: G. fossarum s.l., areas B, C, 
J-L, stations 9-14; i, k-m, o-s) this method 
yielded an estimate that was in contrast with the 
Nm 0 estimate; in 4 cases (Table XI: G. fossarum 
5 ./., G. wautieri and G. orinos; Table XII: G. fossa¬ 
rum s.L , area J, stations o-s) a close agreement of 
values of Nm* and Nm e was found. Although both 
Nm e and Nm*-values are generally low within all 
species and most of the species subsets sampled by 
area, subsampling shows that there are possibly 
high levels of gene flow among many populations. 
In Table XII this is illustrated by G. wautieri (areas 
B and C, stations 12-13, 15) and G. fossarum s.l. 
(area J, stations o-s and p-s) yielding Am^-values 



Fig. 5. CAF gene flow profiles (Slatkin, 1981) in subsets of G. 
fossarum s.l. and G. orinos. For legend see Fig.4. 


of 8.08, 1.17 and 1.75 respectively; Table XIII 
shows that pairwise combinations of individual 
populations may yield high Nm e -\ alues up to 
11.38. Geographic distances in these combinations 
range from 4 to 65 km and in combinations for 
which Nm < 1 from 3 to 55 km. These data indicate 
that even if the possibility of gene flow seems to 
exist between individual populations and areas, 
substantial barriers to genetic exchange may be 
present over small geographic ranges. 

CAF gene flow profiles (Slatkin, 1981; see also 
Caccone, 1985; Caccone & Sbordoni, 1987; Wap- 
les, 1987) of G. p. pulex , G. orinos, G. cf. fossa¬ 
rum , G. wautieri and G. fossarum s.l. are shown in 
Figs. 4, 5. The diagram of Fig. 6 represents typical 
patterns of gene flow (Slatkin, 1981; Waples, 
1987). With reference to this diagram, in Fig. 4, low 
gene flow levels occur in G. orinos, G. cf. fossa¬ 
rum, G. fossarum s.l. 9 and G. wautieri ; and a me¬ 
dium gene flow level in G. p. pulex. All species but 
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Fig. 6. Typical CAF gene flow profiles in species with the indi¬ 
cated levels of gene flow (after Slatkin, 1981; redrawn from 
Waples, 1987). 

G. p. pulex are characterized by high initial p (/)- 
values, and sometimes by wide fluctuations, typical 
for low gene flow species (Caccone & Sbordoni, 
1987). Subsampling (Fig. 5) results in a somewhat 
different pattern. A low to intermediate gene flow 
level is found in G. fossarum s.L from areas B and 
C (stations 9-14); intermediate gene flow levels 
occur in the combined G. fossarum s.L subset from 
areas B, C, J, and L (stations 9—14, i, k-m, o-s) 
and G. orinos , whereas gene flow levels in G. fossa¬ 
rum s.L from area J (stations o-s) and G. fossarum 
s.L from areas K and L (stations i, k-m) are high. 
In these curves, fluctuations are less pronounced. 


Cross-breeding experiments 

The results of the cross-breeding experiments are 
listed in Table XIV. All crosses among populations 
of G. fossarum s.L proved to be completely inter- 
fertile. All crosses involving females of G. fossa¬ 
rum s.L from stations r, 1 and k and males of G. cf. 
fossarum from station b gave FI offspring; all these 
crosses but one (f k x cf. f b ; females of G. fossa¬ 
rum s.L from station k and males G. cf. fossarum 
from station b) failed to produce F2 offspring. In 
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reciprocal crosses, one cross (cf. f b x f r \ females 
G. cf. fossarum from station b and males of G. fos¬ 
sarum s.L from station r) gave FI offspring (one fe¬ 
male), but no F2 was obtained. As for crosses be¬ 
tween populations of G. orinos and G. fossarum 
s.L (including reciprocal crosses), two crosses (o a 
x / 7 , females G. orinos from station a and males 
G. fossarum s.L from station 1; and f k x o a , fe¬ 
males G. fossarum from station k and males G. ori¬ 
nos from station a) gave FI offspring in very low 
numbers; only once (f k x o fl ,), F2 was obtained. 
In crosses between populations of G. orinos and G. 
cf. fossarum , one female from cross cf. f b x o a 
gave FI offspring, but no F2 was obtained. 

Discussion and conclusions 

Genetic differentiation 

1. G. p. pulex & G. p. gallicus 
Inferring genetic relationships from protein poly¬ 
morphisms in the G. pulex- group resulted in some 
unexpected patterns. The topology of genetic struc¬ 
tures was generally determined by geographic dis¬ 
tribution, but sometimes highly incongruent with 
patterns of morphology. Morphologically defined 
taxa covering a large geographic area may show 
little differentiation, while others are highly struc¬ 
tured over much shorter distances. 

Considering the various levels of genetic differ¬ 
entiation among the populations and species 
studied in both UPGMA dendrogram and Wagner 
tree (Figs. 2, 3), the low levels of divergence among 
populations of G. p. pulex and G. p. gallicus are 
striking. Intra-specific /-values of G. pulex (i.e. G. 
p. pulex -l- G. p. gallicus) range from 0.89-0.99 

('mean = ± O' 94 ! Anin ~ Anax = 0. 01-0.11; 

= ±0.06). With reference to other taxa (e.g. 
Scheepmaker et al., 1988; Scheepmaker & Van 
Dalfsen, 1989), /-values of comparisons of geo¬ 
graphically distant populations of G. p. pulex and 
G. p. gallicus may be remarkably high. For in¬ 
stance, G. p. gallicus from station 23 in the 
Pyrenees (southern France; Table I, Fig. 1A, G), 
and G. p. pulex from station t (the Netherlands; 
Table I, Fig. 1A) yield an /-value of 0.93 (D = 
0.07). G. p. pulex has an Eurasiatic distribution 
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(Pinkster, 1972; Karaman & Pinkster, 1977); in 
western Europe it is found from southern Scan¬ 
dinavia to southern France. According to Thiene- 
mann (1950), the extent of the distribution area 
makes a postglacial dispersion highly improbable. 
Therefore, G. p. pulex is to be considered as a 
preglacial, original inhabitant of the lower reaches 
of rivers and streams in the central European 
lowlands as defined by Thienemann (1950). The 
level of divergence between G. p. pulex and G. p. 
gallicus is likely to be of preglacial origin (Figs. 2 
and 3). During the Pleistocene glaciations G. p . gal¬ 
licus (at present restricted to southern France; see 
Table I) probably subsisted in the Mediterranean 
refugial zone, whereas G. p. pulex re-invaded west¬ 
ern European lowlands in postglacial times. 

2. G. wautieri 

Roux (1971) demonstrated that G. wautieri was 
reproductively isolated from G. p. pulex and G. 
fossarum s.L Due to an important overlap of dis¬ 
criminating morphological characters, he had 
difficulties to distinguish the younger individuals 
(especially females) of G. wautieri from G. p. pulex 
and G. fossarum s.L Attempts to support the 
specific status of G. wautieri by other means, such 
as caryological and chromatographical techniques, 
were unsuccessful. Based upon protein polymor¬ 
phisms however, G. wautieri forms a well-defined 
cluster, clearly distinct from G. p. pulex and G. fos¬ 
sarum s.L (Figs. 2, 3). In the present study, only G. 
wautieri from the French Alps was involved. As G. 
wautieri is also known from other areas (e.g. the 
Massif Central in central France, and N.E. of the 
Alps) the data presented here are incomplete and 
may depict a biased intraspecific genetic variation 
pattern. 

5. G. fossarum s.L 

Morphologically very similar (and generally con-' 
sidered conspecific) forms of G. fossarum s.L , G. 
orinos and G. cf. fossarum are occasionally group¬ 
ed separately in distinct clusters (UPGMA and 
Wagner dendrogram; Figs. 2, 3). This is in contrast 
with G. p. pulex , G. p. gallicus , and G. wautieri . 
Subdivision of morphologically very similar forms 
belonging to G. fossarum s.L was discussed by 


Scheepmaker & Van Dalfsen (1989). The material 
studied by these authors (Fig. 1, areas J, K, L; 
Table I, stations i, k-m, o-s) has been completed 
by additional sampling in areas B, C, F and I. In the 
combined data set G. fossarum s.L from areas J, K, 
L (except for station j), B/C, I, and stations 17 and 
o (Table I; Fig. 1A) forms a highly differentiated, 
but monophyletic group (Figs. 2, 3). 

Because of its present day distribution, G. fossa¬ 
rum s.L is considered a preglacial inhabitant of the 
Alps (a “southern glacier border form” according 
to Thienemann, 1950). It is likely that it reached the 
periphery of its distribution as recorded in this study 
(e.g. areas J and K in Fig. 1 A) only since the glacial 
ages. This theory is supported by the high levels of 
genetic similarity among G. fossarum s.L popula¬ 
tions from areas J, C, L and K (Fig. 1 A, Figs. 2, 3). 

4. G. orinos & G. cf. fossarum 
G. orinos and G. cf. fossarum are morphologically 
closely related forms (Scheepmaker et al., 1988). 
The subdivision of both G. orinos (areas D, E, sta¬ 
tion 3-6, vs. area F, stations a, 1-2) and G. cf. fos¬ 
sarum (area H, stations b-d, vs. area G, stations 7, 
8) coincides with an Atlantic-Mediterranean 
watershed (Figs. 2, 3). This subdivision results in 4 
genetically more or less equidistant population sets. 
Levels of genetic differentiation among these sub¬ 
sets (/ = 0.63-0.87; D = 0.14-0.39) are about 
the same as between each of these subsets and sub¬ 
sets of G. fossarum s.L or G. stupendus (I = 
0.73-0.83 and 0.74-0.84, respectively). Cross¬ 
breeding experiments show that G. cf. fossarum 
from area H (station b) and G. orinos from area F 
(station a) are reproductively isolated (Table XIV; 
I = 0.67-0.87). Unfortunately, G. orinos from 
areas D/E and G. cf. fossarum from area G were 
not involved in these crosses. However, crosses of 
G. orinos (station a, area F) and G. cf. fossarum 
(area H, station b) with G. stupendus (I = 0.82 and 
0.75, respectively) and G. fossarum s.L (stations k, 
1, r; / = 0.81, 0.80, 0.84 and 0.69, 0.70, 0.70, 
respectively) were not very successful (Table XIV). 
Differences between these taxa and G. orinos from 
areas D/E, and G. cf. fossarum from area G are in 
the same order of magnitude. G. orinos from areas 
D/E and F (yielding the highest /-values recorded 
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Table XIV. Results of cross-breeding experiments between samples of G. fossarum s.L, G. orinos and G. cf. fossarum. f r = G. fossa- 
rum s.l. from station r;f { = idem, from station \\f k = idem, from station k; o Q = G. ormos from station a; cf./ 6 = G. ct. fossarum 
from station b; legend stations cf. Table I. *)= crosses by Pinkster & Scholl (1984). 


Cross 

9 9 X or or 

Number of 9 9 

Number of cr O’ 

Number of 
ovigerous 9 9 

Number of ovig. 

9 9 with offspring 

F2 

obtained 

fr 

*fr 

32 

32 

15 

10 

+ 

ft 

X// 

126 

122 

81 

18 

+ 

fk 

Xfk 

102 

107 

51 

17 

+ 

o a 

x o a 

107 

77 

60 

12 

+ 

cf./* 

X Cf./* 

167 

114 

93 

20 

+ 

/, 

X fr 

200 

135 

67 

23 

+ 

ft 

X/* 

114 

115 

65 

26 

+ 

ft 

X o a 

57 

47 

38 

- 

- 

ft 

X cf. /* 

163 

112 

89 

1 

- 

fk 

X fr 

192 

119 

117 

27 

+ 

fk 

x cf./* 

114 

89 

69 

11 

+ 

fk 

X o a 

122 

99 

83 

2 

+ 

fk 

O a 

Xf, 

Xf*) 

61 

61 

51 

18 

+ 

O a 

xf, 

82 

52 

12 

1 

- 

°a 

xf k 

104 

104 

45 

- 

- 

°a 

X cf. fu 

197 

214 

131 

1 

- 

cf./* xf r 

230 

190 

34 

1 

- 

cf .ft xf, 

210 

182 

83 

- 

- 

cf-/* X f k 

155 

126 

99 

- 

- 

cf./* 

x o a 

199 

144 

72 

- 

- 

fr 

X// 

142 

87 

110 

37 

+ 

fr 

fr 

X /* 

X 0/) 

248 

100 

182 

46 

+ 

fr 

X cf ,f b 

368 

207 

249 

12 

- 


among subsets of G . orinos) and G. cf. fossarum (I 
= 0.67-0.87), possibly evolve as independent 
units, but they are probably still to be considered 
conspecific. On the other hand, G. cf. fossarum 
from areas G and H (/ = 0.63-0.69) are likely to 
represent distinct taxa, apart from G. orinos , G. 
fossarum s.l . and G. stupendus. However, subsets 
of G. gauthieri from the Iberian peninsula yielding 
/-values as low as 0.56-0.65 were shown to be at 
least partially interfertile (and thus possibly con- 
specific; Scheepmaker et al., 1988). 

5. G. stupendus 

Subdivision in G. stupendus was discussed in detail 
by Scheepmaker (1987). The three population sam¬ 
ples analyzed here (stations 19, 20 and e, belonging 
to G. stupendus forma A, B and C respectively, ac¬ 
cording to Pinkster, 1983) represent morphologi¬ 
cally different and genetically rather discontinuous 
forms (Scheepmaker, 1987). Morphologically, G. 


stupendus forma A is very similar to G. fossarum 
s.l. ; G. stupendus forma B resembles G. cf.fossa- 
rum, and G. stupendus forma C resembles G. ori¬ 
nos in many respects (Table III). 

6. G. fossarum s.str., G. orinos, G. cf. fossarum, 
G. stupendus , G. wautieri vs. G. /?. pulex & G. /as- 
sarum s.L 

G. and G. fossarum s.L are both con¬ 

sidered to be preglacial inhabitants of western 
Europe. Consequently, their level of divergence 
(Fig. 2) must be of preglacial origin. As G. fossa¬ 
rum s.L is considered as a preglacial inhabitant of 
the Alps, inter-areal differentiation (e.g. areas B 
and C vs. area J; Fig. 2) must have occurred since 
the Pleistocene glaciations. From these data we 
may conclude, that species and forms as G. fossa¬ 
rum s.L, G. fossarum s.str., G. orinos, G. cf. fossa¬ 
rum, G. stupendus and G. wautieri, having an 
“earlier” level of divergence, are all likely to be of 
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preglacial origin. These levels do not necessarily 
correspond to those in the Wagner network of Fig. 
3, which is only arbitrarily rooted by the midpoint 
rooting procedure (Farris, 1972); in fact, as will be 
discussed in the next section, there is evidence for 
this network to be rooted at a different inter- 
nodium. 

Unlike G. fossarum 5./., which is widely dis¬ 
tributed in central and northern France, but con¬ 
fined to the area covered with tundras (Brinkman, 
1959) during the last glacial period, G. orinos , G. 
cf. fossarum and G. stupendus have limited dis¬ 
tribution areas in the old Hercynian massifs in 
southern France. According to Figs. 2 and 3, 
smaller geographical units (e.g. G. orinos from 
areas D and E, G. cf .fossarum from area G and H) 
partition a greater portion of the total genetic varia¬ 
tion than the larger ones (e.g. G. p. pulex, G. fossa¬ 
rum s.l. from areas B, C, J, K and L). The smaller 
units correspond to the former Mediterranean refu- 
gia during the Pleistocene glaciations. In contrast 
to the tundra belt, environmental “lowland” con¬ 
ditions may still have prevailed in lower reaches of 
rivers and tributaries. G. orinos 9 G. cf. fossarum 
and G. stupendus generally inhabit the upper 
reaches of rivers. As upper reaches are generally 
separated from each other by lower reaches, this 
may have prevented them from a wider dispersal. 

Phylogenetic relationships among the species and 
forms studied 

The topologies of the dendrograms obtained by the 
UPGMA method and the distance Wagner proce¬ 
dure (Figs. 2 and 3) are not congruent. They con¬ 
trast primarily by suggesting different relationships 
among species and forms; however, the main 
clusters they distinguish agree quite well. 

The UPGMA method and Nei’s (1972) distance 
coefficient require both the assumption of homo¬ 
geneous evolutionary rates. Moreover, Nei’s dis¬ 
tance is unsuitable for deriving phylogenies because 
of its non-metricity (e.g. Farris, 1981). As the dis¬ 
tance Wagner method (Farris, 1972) does not make 
any assumption about evolutionary rates, this 
method seems more appropriate to infer phylo¬ 
genetic relationships. However, the midpoint root¬ 


ing (Farris, 1972) as applied in Fig. 3, also assumes 
homogeneous evolutionary rates. The complete 
data set of the western European members of the G. 
/n/fex-group (i.e. the data from this study linked to 
data sets from the Iberian peninsula and Morrocco; 
Scheepmaker et al., 1988; Scheepmaker, 1990) in 
fact suggests that the Wagner network should not 
be rooted at the midpoint, but at the internodium 
between G. fossarum s.str. and G. fossarum s.l. 
from station j (Fig. 3). In a forthcoming paper, the 
phylogeny of western European members of the G. 
pulex-g roup, derived from protein polymorphisms, 
will be discussed in more detail. 

Comparison of methods used for analyzing popula¬ 
tion structure and estimating levels of gene flow 

The methods used to assess levels of gene flow 
differ in the assumption they depend on. The use of 
F-statistics assumes that genetic divergence by 
stochastic processes is neutralized by the effects of 
gene flow, and not by other processes, such as 
natural selection. Slatkin’s methods (1981, 1985) do 
not depend on the assumption of selective neutrali¬ 
ty (Waples, 1987). Unlike F S7 -values according to 
Wright (1931, 1951, 1965), the 0-values of Weir & 
Cockerham (1984) do not make assumptions con¬ 
cerning numbers of populations, sample sizes or 
heterozygote frequencies. Therefore, according to 
these authors, they are suited to small data sets. 
A r mf - and Am^-values are closely correlated, 
yielding a correlation coefficient (r) of 0.99, 
Nmp sT overestimates levels of gene flow with 
regard to Nm e , possibly due to incorrect assump¬ 
tions about sample sizes, etc. (however, in simula¬ 
tions of Slat kin & Barton, 1989, both methods per¬ 
formed equally well for moderate levels of gene 
flow). 

In all but one cases there is a rather close agree¬ 
ment between Nm*- and Afo^-values. When the 
deviating A7tt*-value is omitted, Nm* and Nm $ 
yield a correlation coefficient of 0.97; for estimates 
that were not corrected for sample size, this value 
is 0.93. The deviating high Nm *-value (G. fossarum 
s.L ; Table XII) can be attributed to the greater 
number of samples involved (about three times 
more than in the other subsets). In fact, the in- 
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fluence of high frequencies of p (1) in areas B and 
C are neutralized by averaging them with the much 
lower p (l)-frequencies in areas J, K and L; 
moreover, as the number of populations increases, 
a number of alleles that were private may not be pri¬ 
vate anymore. One might argue that the complete 
data set of G.fossarum s.l . (19 populations; G.fos - 
sarum s. str is not included) yields a very low 
Nm*- value; however, this value is strongly affected 
by some p (l)-values of 1, due to genetic isolates. It 
may be concluded that we should only compare 
Aw* estimates of data sets of about the same num¬ 
ber of populations. 

As stressed by Caccone & Sbordoni (1987), Slat- 
kin’s (1985) method is appropriate for identifying 
cases in which the majority of populations sampled 
are exchanging (or recently exchanged) genes at a 
rather high rate, but a few populations are genetic 
isolates. 

In many sets and subsets of population samples, 
the 0-value of many polymorphic loci is significant¬ 
ly greater than 0 (Tables VIII-X). It is unlikely, 
that all those loci are affected by natural selection 
at the same time. The correlation between Nm e and 
Am * (presumed not to be influenced by natural 
selection; Barton & Slatkin, 1987) may be in favour 
of this view. As pointed out by Slatkin (1985), a 
basic problem (which, although for different rea¬ 
sons, also affects Nm estimates obtained with F - 
statistics) is the effect on the estimated value of 
Nm * when migration is between nearest neighbours 
in a stepping stone model, and the populations sam¬ 
pled are separated by very large distances (e.g. G. 
p. pulex). In that case, an allele that is encountered 
in only one population could be present in several 
adjacent populations that were not sampled. A 
minimum number of samples, relative to the area 
covered and its structure and to the presumed dis¬ 
persal abilities of the species under study, is recom¬ 
mended. Spatial heterogeneity may also be a source 
of bias in estimating levels of gene flow. Splitting 
up the dataset of a species in different subsets also 
contributes to a more accurate Nm estimate. 
Another problem is the occurrence of bottlenecks. 
When gene flow levels are low (Nm < 1), rare al¬ 
leles will disappear due to genetic drift and bottle¬ 
necks. Consequently, Nm * will overestimate the ac- 
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tual levels of gene flow, representing rather histori¬ 
cal patterns of genetic exchange. 


Population structures of morphological species 

G. p . pulex and G. p . gallicus show little genetic 
differentiation, although the geographical area 
covered by these forms is large. The gene flow pro¬ 
file of G. p. pulex (Fig. 4) resembles the characteris¬ 
tic curve of high gene flow species. However, G. p. 
pulex is significantly subdivided at most loci (Table 
IX) and estimated gene flow levels (Nm* and Nm e ) 
are low (less than one; Table XI). The reasons may 
be twofold. Firstly (probably the most important), 
due to the small number of samples (and small sam¬ 
ple sizes) relative to the large area covered, a basi¬ 
cally clinal variation pattern will appear discontinu¬ 
ous, and consequently, estimated levels of gene 
flow will be low. Secondly, low levels of genetic 
differentiation and low gene flow profiles (CAF) 
may be more indicative for historical patterns of 
genetic exchange than of current patterns, whereas 
the present distribution of allelic variants indicates 
that genetic exchange in the past may have been 
higher. 

Contrary to G. p. pulex , populations of G. fos - 
sarum s.l., G. orinos , G. cf. fossarum and G. 
stupendus in fact are significantly subdivided at 
many polymorphic loci. Significant subdivision 
even occurred at higher levels than those recorded 
in cave-dwelling species (Caccone & Sbordoni, 
1987), and estimated levels of gene flow (Aw* and 
Aw 0 ) are generally smaller than one. An explana¬ 
tion is possibly provided by the (in)stability of the 
environment. Although isolated, the cave environ¬ 
ment (and consequently population sizes) may be 
stable over many years. On the other hand, streams 
and rivers (especially in the Alps and the Mediterra¬ 
nean area) may be subject to dramatic environmen¬ 
tal changes, such as almost complete desiccation 
(Scheepmaker, 1987) or strongly increased dis¬ 
charge and subsequent erosion due to heavy rain¬ 
fall. Repeatedly, surviving population fragments 
are subjected to stochastic processes. This hypothe¬ 
sis is supported by area dependent gene flow levels 
in G.fossarum s.l. (Table XII). Levels of gene flow 
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in the French Alps (areas B, C) are low; in northern 
France however (area J), where moderate environ¬ 
mental conditions prevail, levels of gene flow are 
higher (Nm* and Nm e are substantially > 1). As a 
matter of fact, Goedmakers (1981) and Goedma- 
kers & Pinkster (1981), who demonstrated that 
gammarids are actively and sometimes massively 
migrating organisms, studied G. fossarum s.L in 
this very area. However, this activity may be quite 
different from one drainage system to another 
(Pinkster, pers. comm.). The observed population 
subdivision in structured environments agrees with 
findings of Gooch & Hetrick (1979) and Scheep¬ 
maker (1987), who studied freshwater gammarids 
in karst areas; these authors also reported sharp dis¬ 
continuities in population structure. 


Estimated levels of gene flow and environmental 
structure 

The results reported here confirm that gene flow oc¬ 
curs primarily in relation to environmental struc¬ 
ture. For example, in the highly structured French 
Alps (areas B, C) the Nm * estimate within G. fossa- 
rum s.L is 0.20; in area J, a more uniform area, this 
value is 1.51. 

For comparison, within the genus Gammarus , 
between two population samples of the freshwater 
form G. duebeni celticus in Brittany (France) the 
Nm e estimate is 0.60; between two population sam¬ 
ples of the euryhaline brackish form G. d . duebeni 
from area J (Fig. 1A) and Petten (North-Holland) 
this estimate yields 5.98 (Scheepmaker, 1990); in a 
marine environment, Nm e estimates of G. salinus 
(18 samples, ranging from Denmark to southern 
Brittany) and subsets of G. zaddachi (7 and 11 sam¬ 
ples, ranging from Tromso to Brittany) calculated 
from data provided by Bulnheim & Scholl (1981) 
yield 4.26, 7.79 and 1.45, respectively; in a lowland 
and “polder’’ environment, G. tigrinus (10samples 
from Holland to Denmark; calculated from Buln¬ 
heim, 1985) yields an Nm e estimate of 1.56 
(Scheepmaker, unpublished data). These data cor¬ 
roborate the suggested relation between environ¬ 
mental structure and estimated levels of gene flow. 


Genetic differentiation and cross-breeding ex¬ 
periments 

Cross-breeding experiments by Roux (1971), Goed¬ 
makers (1972), Pinkster (1972, 1983), Pinkster & 
Scholl (1984), Scheepmaker (1987), and results 
from crosses presented in this study (Table XIV) 
suggest that clusters formed by G. pulex , G. wau - 
tieri, G. st upend us, G. fossarum s.L from areas J, 
K, L, B, C, F and I (Fig. 1A), G. orinos from area 
F, and G. cf .fossarum from area H, act as indepen¬ 
dently evolving evolutionary units. Forms within 
those clusters are thus to be considered conspecific. 
The position of G. fossarum s.L from stations n 
and j, G. fossarum s. str. y G. orinos from areas D 
and E, and G. cf. fossarum from area G remains 
uncertain. Although not morphologically distinct 
from other forms, it seems quite well possible that 
these forms also evolve independently, ultimately 
resulting in distinct species. 
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